Standard random walks and trapping on the Koch network 
with scale-free behavior and small-world effect 
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A vast variety of real-life networks display the ubiquitous presence of scale-free phenomenon and 
small- world effect, both of which play a significant role in the dynamical processes running on 
networks. Although various dynamical processes have been investigated in scale-free small-world 
networks, analytical research about random walks on such networks is much less. In this paper, 
we will study analytically the scaling of the mean first-passage time (MFPT) for random walks on 
scale-free small-world networks. To this end, we first map the classical Koch fractal to a network, 
called Koch network. According to this proposed mapping, we present an iterative algorithm for 
generating the Koch network, based on which we derive closed-form expressions for the relevant 
topological features, such as degree distribution, clustering coefficient, average path length, and 
degree correlations. The obtained solutions show that the Koch network exhibits scale-free behavior 
and small-world effect. Then, we investigate the standard random walks and trapping issue on the 
Koch network. Through the recurrence relations derived from the structure of the Koch network, we 
obtain the exact scaling for the MFPT. We show that in the infinite network order limit, the MFPT 
grows linearly with the number of all nodes in the network. The obtained analytical results are 
corroborated by direct extensive numerical calculations. In addition, we also determine the scaling 
efficiency exponents characterizing random walks on the Koch network. 

PACS numbers: 05.40.Fb, 89. 75. He, 05.60.Cd, 05.10.-a 



I. INTRODUCTION 

Complex networks have been acknowledged as an in- 
valuable tool for describing real-world systems in nature 
and society P, [3, H, 0, Q . Extensive empirical studies 
have uncovered that a lot of real networks share several 
remarkable features One of the most relevant is the 
scale-free behavior, that is, various real networks exhibit 
a power-law degree distribution P{k) ~ 0- Another 
very important observation is that most real-life systems 
are characterized by ubiquitous small- world effect Q , in- 
cluding large clustering coefficient @ and small average 
path length (APL) [l^]- Both the scale- free behavior 
and the small-world effect have a profound impact on 
almost all dynamical processes taking place on the net- 
works d i, i, i [m, E [n, [11 [n, m [13. 

Among various dynamical processes, random walks on 
networks are fundamental to many branches of science, 
and have received considerable attention from the scien- 
tific community [H [H, HO, [SI, HI HI, [H, H. As a 
primary dynamical process, random walks are related to 
a plethora of other dynamics such as transport in me- 
dia [2&I , disease spreading [271 , and target search [H, |2§| 
to name a few. On the other hand, random walks are use- 
ful for the study of topological structure (e.g. between- 
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ness and average path length [l§ . [25j ) and community 
detection [13] on networks. In particular, as an integral 
theme of random walks, trapping is related to a wide 
variety of contexts 31], such as pho ton- harvesting pro- 
cesses in photosynthetic cells j32l. T33| and characterizing 
similarities between the elements of a database ^] . It is 
thus of theoretical and practical interests to study stan- 
dard random walks and trapping problem on complex 
networks. 

The main interesting quantity closely related to ran- 
dom walks is the mean first-passage time (MFPT), and 
a central issue in the study of random walks is how the 
MFPT scales with the size of the system [H, [H, [11]. 
In the past several years, a lot of endeavors have been 
devoted to studying the intrinsic relations between the 
scaling behavior of random walks and the underlying 
topological structure of networks [3, [l^, [Sll, [sj ■ The re- 
sults of these investigations uncovered many unusual and 
exotic features of complex networks, especially of small- 
world and scale-free networks. In spite of their useful 
insight, most of previous jobs focused on numerical in- 
vestigation, the analytical results of MFPT (in particular 
the average for all pairs of nodes) for standard random 
walks and trapping problem have been far less reported, 
with the exception of some graphs with simple topology, 
such as regular lattices (s^l, Sierpinski fractals [38, 39], T 
fractal [40[, and deterministic scale- free trees [4l|, [Uf, as 
well as other structures [20I [2ll . [23 ] , and exhaustive an- 
alytical research on scale-free and small-world networks 
with loops is still missing. 

In this paper, we analytically investigate the scaling 
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FIG. 1: (Color online) The first two generations of the con- 
struction for Koch curve. 



behavior of MFPT on complex networks with scale-free 
phenomenon and small- world effect. To achieve this goal, 
we first propose a relevant deterministic network, named 
Koch network, which is based on the classical fractal — 
Koch curve. We then suggest a minimal iterative algo- 
rithm generating the Koch network, on the basis of which 
we give in detail a scrutiny of the network architecture. 
The analysis results show that the Koch network is si- 
multaneously scale free, small world, and has large clus- 
tering coefficient. Particularly, we show that the Koch 
network is completely two-point uncorrelated by exactly 
computing two relevant quantities of two-point correla- 
tions, which has never been previously found in other de- 
terministically growing networks. Eventually, we study 
the scalings of MFPT on the proposed network for the 
standard random walks, together with a special case of 
trapping problem with the trap fixed at a hub node. We 
present that the MFPT averaged over all couples of nodes 
scales linearly with the network order (total number of 
nodes). We derive a closed- form solution for the MFPT 
characterizing the trapping process, which also grows as 
a linear function of the network order. We also compare 
the our results with those previously obtained for other 
scale-free networks. 



II. NETWORK CONSTRUCTION 
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FIG. 2; (Color online) Construction for the Koch network 
and labels of its nodes, showing the first two steps. 

the segment upon which the new triangle is established. 
Thus, we get Si. For each line segment in 5*1, trisect it 
and draw an equilateral triangle based on the resultant 
middle small segment to obtain 5*2. Repeat recursively 
the procedure of trisection of existing line segments in 
last generation and addition of triangles. In the infinite t 
limit, we obtain the famous Koch curve St, whose Haus- 
dorff dimension is df = SJi^ [45|. In Fig. [TJ we show 
schematically the structure of S'2. fn fact, this fractal 
can be easily generalized to other dimensions [4^ . 

From the Koch curve we can easily construct a net- 
work (called Koch network) using a simple mapping as 
follows. In the Koch network, nodes (vertices and sites) 
correspond to the sides (excluding those deleted) of the 
triangles constructed at all generations of the Koch curve 
as shown in Fig. [TJ That is to say, for every triangle cre- 
ated at some generation, its two newly-born sides are 
mapped to two nodes, while the removed side is not. We 
make two nodes connected if the corresponding two sides 
of the Koch curves contact each other. For uniformity, 
the three sides of the initial equilateral triangle of Sq also 
correspond to three different nodes. Note that after the 
birth of each side of a triangle constructed at a given gen- 
eration, although some segments of it will be deleted at 
subsequent steps, we look upon its remaining segments 
as a whole and map it to only one node. Figure [5] shows 
the network associated with 6*2. 



The network under consideration is derived from the 
Koch curve. To define the network, we first introduce the 
classical fractal, Koch island, also known as Koch curve, 
Koch snowflake or Koch triangle, which was proposed by 
Koch [43*1. This well-known fractal denoted by St after 
t generations is constructed as follows [44]. Start with 
an equilateral triangle and denote this initial configura- 
tion as 5*0. Perform a trisection of each side of this ini- 
tial triangle and construct an equilateral triangle on each 
middle segment, so that the interior of the added trian- 
gle lies in the exterior of the base triangle, then remove 



III. GENERATION ALGORITHM OF KOCH 
NETWORK 

According to the construction process of Koch curve 
and the proposed mapping from Koch curve to Koch net- 
work, we can introduce with ease an iterative algorithm 
to create Koch network, denoted by K{t) after t gener- 
ation evolutions. The algorithm is as follows. Initially 
(t = 0), K{Q) consists of three nodes forming a trian- 
gle. Then, each of the three nodes of the initial triangle 
gives birth to two nodes. These two new nodes and its 
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FIG. 3; (Color online) Iterative construction method for the 
network. 



mother node are linked to each other shaping a new tri- 
angle. Thus we get K{1) (see Fig. E]). For t > 1, K{t) 
is obtained from K{t — 1). We replace each of the exist- 
ing triangles of K{t — 1) with the connected cluster on 
the right-hand side (rhs) of Fig. [3] to obtain K{t). The 
growing process is repeated until the network reaches a 
desired order. Figure[2]shows the network growth process 
for the first two steps. 

Next we compute the order and the size (number of 
all edges) of Koch network K{t). To this end, we first 
calculate the total number of triangles existing at step t, 
which we denote as L^(t). By construction, this quan- 
tity increases by a factor of 4, i.e., L^it) = 4LA(t — 1). 
Considering the initial condition £a(0) = 1, it follows 
that LA{t) = 4*. Let Lv{t) and Le{t) be the respective 
number of nodes and edges created at step t. Notice that 
each triangle in K{t — 1) will lead to an addition of six 
new nodes and nine new edges at step t, then one can eas- 
ily obtain the following relations: L„(t) = 6La(^ ^ 1) = 
6 X 4*-i and Le{t) = ^Li^{t - 1) = 9 x 4*^^ for arbitrary 
t > Q. From these results, we can compute the order and 
the size of Koch network. The total number of vertices 
Nt and edges Et present at step t is 



^L„(t,)=2x4* + 1 



(1) 



and 



Et^Y. = 3x4*, 

respectively. Thus, the average degree is 
2Et 6x4* 



(2) 



(fc) 



Nt 2x4* + !' 



(3) 



which is approximately 3 for large t, showing that Koch 
network is sparse as most real systems. 



A. Degree distribution 

We define ki{t) as the degree of a node i at time t. 
When node i is added to the network at step ti {ti > 0), 
it has a degree of 2, viz., ki{ti) — 2. To determine ki{t), 
we first determine the number of triangles involving node 
i at step t that is represented by L^{i, t). These triangles 
will create new nodes connected to the node i at step 
t + 1. Then at step i^, L\{i,ti) — 1. By construction, 
L^iht) — 2Li^(i,t — 1). Since L^{i,ti) — 1, one can 
derive Li\{i,t) = 2*~*'. Note that the relation between 
ki{t) and LA{i,t) satisfies: 



fc,(i) = 2iA(«,t) = 2*+i 



(4) 



In this way, at time t the degree of node i has been com- 
puted explicitly. From Eq. (|3]), one can see that at each 
step the degree of a node doubles, i.e.. 



hit) = 2k,{t-l). 



(5) 



Equation ([4]) shows that the degree spectrum of Koch 
network is discrete. It follows that the cumulative degree 
distribution ^ is given by 



Pcnm{k) = ^ ^L,„(r)=^''^' 

* T<t, 



1 



2 x 4* + 1 



(6) 



Substituting for ti in this expression using ti ~ t+1 — |^ 
gives 

, , 2 X 4* X 4 X fc-(ln4/ln2) _^ i 

P^Uk) . ^^^^ . (7) 

When t is large enough, one can obtain 

Pcnmik) - 4 X fc-2. (8) 

So the degree distribution follows a power-law form with 
the exponent 7 = 3. Note that this exponent of degree 
distribution is the same as that of the Barabasi- Albert 
(BA) model ^7,]- 

Before closing this section, we compute another quan- 
tity (fc^), i.e., the fluctuations of the connectivity distri- 
bution, which is useful for the calculation of Pearson cor- 
relation coefficient that will be discussed in the following 
text. The quantity (fc^) is given by 



(9) 



u=o 



IV. STRUCTURAL PROPERTIES OF KOCH 
NETWORK 

Now we study some relevant characteristics of Koch 
network K{t), focusing on degree distribution, clustering 
coefficient, average path length, and degree correlations. 



where k{ti,t) is the degree of a node at step t which 
was generated at step ti. Combining previously obtained 
results, we find 



2 2 X 4*(3t-|-6) 



2 X 4* + 1 



(10) 
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FIG. 4: Semilogarithmic plot of average clustering coefficient FIG. 5: Average path length dt versus network order Nt on a 
Ct versus network order Nt. semilogarithmic scale. The solid line is a guide to the eye. 



B. Clustering coefficient 

The clustering coefRcient [1| of a node i with a degree 
ki is given by Ci = 2ei/[ki{ki — l)], where is the number 
of existing triangles attached to node i, and ki{ki — l)/2is 
the total number of possible triangles including i. Using 
the connection rules, it is straightforward to calculate an- 
alytically the clustering coefficient C(fc) for a single node 
with degree k. In the preceding section, we have obtained 
2ei = ki for all nodes at all steps. So there is a one-to- 
one correspondence between the clustering coefficient of 
a node and its degree. For a node of degree k, we have 



C{k) = 



1 



(11) 



which is inversely proportional to k in the limit of large 
k. The scaling of C{k) ~ k~^ has been observed in many 
real- world scale- free networks [46^1 . 

After t generation evolutions, the clustering coefRcient 
Ct of the whole network, defined as the average of Ci s 
over all nodes in the network, is given by 



Ct 



1 * 

— T 

Nt ^ 



1 



Mr - 1 



(12) 



where the sum runs over all the nodes and M,. is the 
degree of those nodes created at step r, which is given 
by Eq. (|4|). In the limit of large Nt, Eq. (fT2)) converges 
to a nonzero value C = 0.82008, as shown in Fig. U) 
Therefore, the Koch network is highly clustered. 



of K(t) as a matrix in which the entry dij is the shortest 
distance from node i to node j, then dt is defined as the 
mean of dij over all couples of nodes, 



dt = 



Dt 



Nt{Nt - l)/2 ' 



where 



Dt = 



E 



(13) 



(14) 



i<iK(t), jeK(t) 



denotes the sum of the shortest path length between two 
nodes over all pairs. It should be mentioned that in 
Eq. (fTi]) . for a couple of nodes i and j (i ^ j), we only 
count dij or dji, not both. In the Appendix, we pro- 
vide the detailed derivation for the APL. The obtained 
analytical expression for dt is 



4+ 14 X 4* + 12t X 4* 
3 (4* -t- 1) 



(15) 



which approximates At in the infinite t, implying that the 
APL shows a logarithmic scaling with the network order. 
Therefore, the Koch network exhibits a small-world be- 
havior. We have checked our analytical result against 
numerical calculations for different network orders up to 
t = 10 which corresponds to A^io — 1 048 577. In all the 
cases we obtain a complete agreement between our theo- 
retical formula and the results of numerical investigation 
(see Fig. [SJ. 



C. Average path length 

Let dt denote the APL of the Koch network K{t). 
Since the Koch network is self-similar, the APL can be 
computed analytically to obtain an explicit formula, by 
using a method similar to but different from those in 
Refs. [13,14^. We represent all the shortest path lengths 



D. Degree correlations 

Degree correlations are a particularly interesting sub- 
ject in the field of network science [i^, because they can 
give rise to some interesting network structure effects. An 
interesting quantity related to degree correlations is the 
average degree of the nearest neighbors for nodes with 
degree fc, denoted as knn{k), which is a function of node 
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degree k [50|- When fcnn(fc) increases with k, it means 
that nodes have a tendency to connect to nodes with 
a similar or larger degree. In this case the network is 
defined as assortative [51]. In contrast, if fcnn(fc) is de- 
creasing with fc, which imphes that nodes of large degree 
are likely to have near neighbors with small degree, then 
the network is said to be disassortative. If correlations 
are absent, fcnn(fc) = const. 

We can exactly calculate k^nik) for Koch network us- 
ing Eqs. ([U and ([5]) to work out how many links are made 
at a particular step to nodes with a particular degree. By 
construction, we have the following expression [s^ . [53| : 



^nn {k) 



L4U)k{U,t) 



t'=t 



' t'-=ti-l 

J2 L,{t'Mt[,U-l)k{t[,t) 



t'=o 



L,{U)k{U,t[-l)k{t'M +1 (16) 



t'^=ti+i / 

for k — 2*+i~*'. Here the first sum on the right-hand 
side accounts for the links made to nodes with larger 
degree (i.e., t[ < U) when the node was generated at ti. 
The second sum describes the links made to the current 
smallest degree nodes at each step > ti. The last term 
1 accounts for the link connected to the simultaneously 
emerging node. After some algebraic manipulations, we 
obtain exactly 



(17) 



Therefore, two node correlations do not depend on the 
degree. On the other hand, Eq. P7| shows that for 
large t, k-a^{k) is approximately a logarithmic function 
of the network order Nt, namely knn{k) ~ In iVt. Note 
that the same behavior has also been observed in the BA 
model nil. 

Degree correlations can be also described by a Pearson 
correlation coefficient r of degrees at either end of a link. 
It is defined as 0, [sl, HI 



(fc)(fc^fcnn(fc)) - {ky 
{k){k^) - (fc2)2 



(18) 



If the network is uncorrelated, the correlation coefRcient 
equals zero. Disassortative networks have r < 0, while 
assortative graphs have a value of r > 0. Substituting 
Eqs. (O, dTO]), and ^ into Eq. HH), we can easily see 
that for arbitrary t > 0, the numerator of Eq. (|18p is 
always equal to zero. Thereby, r also equals zero, which 
again indicates that Koch network shows the absence of 
degree correlations. 



V. RANDOM WALKS 



As addressed in Sec. lIVl the Koch network exhibits ex- 
clusive topological properties not simultaneously shared 



by other networks. Thus, it is worthwhile to study dy- 
namical processes occurring on the network. In this 
section we consider simple random walks on the Koch 
network defined by a walker such that at each step the 
walker, located on a given node, moves to any of its near- 
est neighbors with equal probabilities. 



A. Scaling efficiency 



We follow the concept of scaling efficiency introduced 
in [4l|. Denote by Tij the first-passage time (FPT) be- 
tween two nodes i and j in a network. Let Tu be the 
mean time for a walker returning to a node i for the first 
time after the walker has left it. When the network or- 
der grows from N to gN ^ one expects that in the infinite 
hmit of N 



T,,igN)^g''^T,,{N), 



(19) 



where Sij is defined as the scaling efficiency exponent. 
An analogous relation for Tu defines an exponent Su. 

One can confine the scaling efficiency in the nodes al- 
ready existing in the network before growth. Let Tlj{gN) 
be the mean first-passage time in the network under con- 
sideration, averaged over the original class of nodes (be- 
fore growth). Then the restricted scaling efficiency expo- 



nent \ij is defined by the relation 

T^igN) ^ g^^^T,,{N). 



(20) 



Similarly, we can define Xu. 

After introducing the concepts, in the following we will 
investigate random walks on the Koch network following 
a similar but obviously different method used in [U \^ . 



B. First- passage time for old nodes 

Consider an arbitrary node i in the Koch network K{t) 
after t generation evolution. Note that for the sake of 
simplicity, we also denote K{t) by Kt, and both deno- 
tations will be used alternatively in the following text. 
From Eq. ([5]), we know that upon growth of the network 
to generation t + 1, the degree ki of node i doubles, that is 
to say, it increases from ki to 2ki. Let the FPT for going 
from node i to any of the ki old neighbors be T and let 
the FPT for going from any of the ki new neighbors to 
one of the ki old neighbors be A. Then we can establish 
the following equations (see Fig. ^ : 



T = 



i + i(l + A), 



A = i(l+T) + i(l + A), 



(21) 



which leads to T = 4. Therefore, the passage time from 
any node i (i G Kt) to any node j {j G ^^t+i) increases 
four times, on average, upon growth of the network to 
generation t + 1, i.e., 



TUN, 



t+i) 



4ry(iVt). 



(22) 
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FIG. 6: (Color online) Growth of trapping time in going from 
Kt to Kt+i- Node i £ Kt has ki neighbor nodes in generation 
t (•) and ki new neighbor nodes in generation t+1 (■). A new 
neighbor of node i has a degree of 2, and is simultaneously 
linked to another new neighbor of i. 

For explanation, see Refs. [l^, |4l| and related references 
therein. Since the network order approximately grows by 
four times in the large t limit [see Eq. ([T|)]. This indicates 
that the scaling efficiency exponent for old nodes is Ay = 
1, which is the same as that of the recursive scale- free tree 
addressed in Ref. [ilj . 

Next we continue to consider the return of FPT to 
node i. Denote by T/j the FPT for returning to node i in 
Kt^i. Denote by Tj^ the FPT from j — an old neighbor 
of i (i e Kt) — to i, in Kt+i- Analogously, denote by Ta 
the FPT for returning to i in Kt and denote by Tji the 
FPT from the same neighbor j, to i, in Kt. For Kt, we 
have 

where fli{t) is the set of neighbors of node i, which be- 
longs to Kt. On the other hand, for Kt+i, 

The first term on the rhs of Eq. (|24p accounts for the 
process in which the walker moves from node i to its new 
neighbors and back. Since among all neighbors of node 
j, half of them are new, which is obvious from Eq. ([S]), 
such a process occurs with a probability of ^ and takes 
three time steps. The second term on the rhs interprets 
the process where the walker steps from i to one of the 
old neighbors j previously existing in Kt and back; this 
process happens with the complimentary probability (1 — 

Using Eq. ([22|l to simphfy Eq. ([24)l . we can obtain 

Tl,^2Tu = A'/^Tu. (25) 

In other words, 

Tl,{Nt+i)=4'^'Tu{Nt). (26) 

Thus, the scaling efficiency exponent Xu = ^, which is 
less than 1. Recall that for the recursive scale- free tree, 



its scaling efficiency exponent is Xu = 1 — j^l < 1 jilj . 
which means that in the Koch network it is more difficult 
for the walker to return to the origin than in the recursive 
scale-free tree, when the networks grow in size. 

C. First-passage time for all nodes 

We continue to compute Tjiji, which is the FPT to 
return to a new node j' G Kt that is a neighbor of node 
i G Kt-i. Notice that when j' was generated, another 
node j" emerged simultaneously, connected to i and j' 
(see Fig. E]). Denote by Ti the FPT from i to j' and 
denote by B the FPT to return to i (starting off from i) 
without ever visiting j' and j" . Then we have 

T,',' = i(l+Ti) + i(l + T,.,-), (27) 



T,.,. = -xl + -(l + Ti), (28) 

and 

Ti = Y + h^ + + + ^i)- (29) 

Equation can be interpreted as follows: with proba- 
bility -J- [ki being the degree of node i in Kt), the walker 
starting from node i would take one time step to go to 
node j'; with probability j-, the walker takes one time 
step to move to node j" then takes time Tjiiji to reach 
j'; and with the remaining probability ^p-^, the walker 
chooses uniformly a neighbor node except j' and j" and 
spends on average time B in returning to i then takes 
time Ti to arrive at node j' . 

In order to close Eqs. (|?7|) and we express the 

FPT to return to i as 

Tu{Nt) - X 3 + X 3 + ^^V^B. (30) 

Eliminating Ti, Tj"j'^ and we obtain 

T,','{Nt)^jTu{Nt). (31) 

Combining Eqs. (gl), and we have 

T,,,,{Nt) = 3xA' ^3Nt. (32) 

Thus, in spite of the fact that simultaneously emerging 
new nodes are linked to different nodes with various de- 
grees, they have the same mean return time. Iterating 
Eqs. (|^ and (PT|) . we have that in Kt there are Ly{e) 
(0 < e < t) nodes with Tu = 2*+l+^ This, together 
with Eqs. ^ and means that for an arbitrary node 
i (born at step ti) with degree ki{t) = 2*+^"*' at time t, 
the FPT to return to i is Tu = Note that a similar 

expression has been obtained for the periodic lattices [s^ 
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and random networks [l8| by using a different derivation 
metfiod. Thus, the mean return time Tn is not reliant 
on the details of the global structure of Koch network. 
It only relies on the network size and the connectivity of 
the node: the larger the degree of the node, the smaller 
the FPT to return. 

Taking the average of Tu over all nodes in Kt leads to 

{Tu)t = ?^^f3x2 



t+i 



2 X 4* + 1 



1 



2 X 4* + 1 V 7 



6x2 



t+2 



6x4x2 



t+2 



+6 X 4*-i x 2*+!+*) 



24 



X 16* 



18 

y 



X 2* 



(33) 



For infinite t, (Tii}t ~ ^ x 4^ ^ Nt, implying that Sa = 
1, an efficiency scaling identical to that of the recursive 
scale- free tree |4ll] . Therefore, the scalings of new nodes 
play a dominant role in the average of mean return time 
for all nodes. 

Next we calculate Ty in Kt, which is FPT from an 
arbitrary node i to another node j. Since each newly 
created node has a degree of 2 and is linked to an old 
node and a simultaneously emerging new node, and these 
three nodes form a triangle, the FPT T^/j from node i' — 
a new neighbor of the old node i — to j equals Ty plus 2 
(i.e., Ti'i) and thus has little effect on the scaling when 
the network order N is very large. Therefore, we need 
only to consider FPT Ty' from i to j' — a new neighbor 
of j , which can be expressed as 



Tiji — Tij 



33' ■ 



(34) 

Suppose that when j' was born, it connected to node j 
and a simultaneously emerging node j" (by construction, 
j" was also linked to node j), then we have 



T , V — 

-^3 3 ^ 



- T, 



and 



(35) 



(36) 



Inserting Eq. ([36|l to Eq. ([35|l . we obtain 



- 2. 



-^33 ~ 2 -^3 3 



Substituting Eq. 
results in 



37]) and Eq. §^ for T. 



T, 



ij' 



It-, 

3 ^ 



'-3'3' 



Nt 



(37) 

into Eq. ^ 



(38) 



where Eq. ([22]) has been used. Therefore, we have 

{T,j)tr^Nt, (39) 

which shows that the mean transit time between arbi- 
trary pairs of nodes is proportional to the network order. 
Equation ((39|) also reveals that 6ij is a constant 1. Notice 
that the linear scaling of the average traverse time with 
network order has been previously obtained by numeri- 
cal simulations for the Apollonian networks [37[ and the 



pseudofractal scale- free we b 141 
wen studied [H, [s^, [H, [Hi, [fi 



both of which have been 
6ii[6l,[6l,[6i[6i,[6i,[63. 



VI. RANDOM WALK WITH A TRAP 

In the preceding section, we have obtained the scaling 
between the MFPT and the network order. Although 
scaling theory is central in studying diffusive process on 
a variety of deterministic media, it has been noted that 
scaling laws do not provide a complete picture of dy- 
namic phenomena on deterministic media, e.g., regular 
fractals [6^, and exact relationships are useful. In this 
section we study the trapping problem of a simple unbi- 
ased Markovian random walk of a particle on network Kt 
in the presence of a trap or a perfect absorber located on 
a given node, which absorbs all particles visiting it. Our 
aim is to give some further insight into the trapping pro- 
cess and the independence of the average trapping time 
(TT) on the size of the underlying system, by providing 
a rigorous solution of the mean trapping time (MTT) as 
a function of the network order. 



A. Formulating the trapping problem 

For the convenience of description, we distinguish dif- 
ferent nodes by labeling all the nodes belonging to Kt in 
the following way. The initial three nodes in Kq are la- 
beled as 1, 2, and 3, respectively. In each new generation, 
only the new nodes created at this generation are labeled, 
while the labels of all old nodes remain unchanged, i.e., 
we label new nodes as M + 1, M + 2, . . ., M + AM, where 
M is the total number of the pre-existing nodes and AM 
is the number of newly created nodes. Eventually, every 
node is labeled by a unique integer, at time t all nodes 
are labeled from 1 to iVt = 2 x 4* -I- 1 (see Fig. [5]). 

Before proceeding further, we give the following def- 
initions. Let Aij be an element of the adjacency ma- 
trix Af of network Kt such that An = 1 if nodes i 

otherwise. 



Lt Oi UeUWOrK J\t ou-^-ii u±icn, js-ij 

and j are connected by an edge and Aij 
Thus, the degree of a vertex i in Kt is dy- = X^j^y- 
And the diagonal degree matrix Dt of Kt is defined by 
T>t = diag(d„i , d^^, ■ ■ ■ , dy., . . . , d^^J. Finally, we define 
Wt = At ■ D^^, where D^^^ is the inverse matrix of Dj, 
then the normalized Laplacian matrix of network Kt is 
lit — It ^ Wt, in which It is an identity matrix with order 
Nt X Nt. 

We locate the trap at node 1 (due to the symmetry, the 
trap can be also located at node 2 or 3, which does not 
have any effect on MFPT), denoted as zt- Note that the 
particular selection we made for the trap location makes 
the analytical computation process (that will be shown 
in detail in the following text) easily iterated as we can 
identify the trap node ix since the first generation. At 
each time step (taken to be unity), the walker selects 
uniformly among the neighbors of the current node (ex- 
cluding the trap) and takes a step to one of them. Since 
node 1 is one of the three nodes with the largest degree, 
it is easily seen that in the presence of the trap ix fixed 
on node 1, the walker will be inevitably absorbed (4l| . 

This trapping process can be described by specifying 
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the set of transition probabilities Wij for the particle of 
going from node i (except the trap z^) to node j. We 
can regard Wij as an element of the matrix Wj , which 
is a submatrix of Wt with the row and the column cor- 
responding to trap being removed. That is to say, 
is a {Nt — l)-order submatrix of Wt with the first row 
and column being deleted. Similarly, one can define , 
, if and LJ , thus we have Wf = AJ ■ (Df )-i and 
Lf = If-Wj. The inverse of L^, {Lf)-\ is the funda- 
mental matrix of the Markovian chain representing the 
unbiased random walk with a trap. 

An interesting quantity related to the trapping process 
is the mean residence time (MKT), which is the mean 
time that a random walker spends at a given site prior to 
being trapped. In fact, the MKT is the mean number of 
visitations of a given site by the walker before trapping 
occurs, and finding the MKT at node i, starting from 
node j, is equivalent to finding the element L~j of the 
fundamental matrix [6^ . 

Another quantity of interest is the TT. In network Kt, 
the trapping time T/ of a given site i is the expected time 
for a walker starting from i to first reach the trap. By 
definition, trapping time Tf is the sum of the MKTs over 
all nodes except ix, i.e.. 

Aft 

t!-J2^^- (40) 

Then, the MTT or the mean first-passage time 
(MFPT), {T)t, which is the average of T/ over all ini- 
tial nodes distributed uniformly over nodes in Kf other 
than the trap, is given by 

, ATt Nt 

i=2 i=2 j=2 

The quantities of TT and MTT are very important 
since they measure the efficiency of the trapping process: 
the smaller the two quantities, the higher the efficiency, 
and vice versa. Equations (|iD|) and ||1T|) show that the 
problem of calculating Tf and {T)t is reduced to finding 
the sum of elements of matrix (LJ)'^. In Tables HI and [Til 
we list separately T/ of some nodes and {T)t for different 
network orders up to t = 6. From Table U one can easily 
observe that for a given node i, the relation T*+^ ^ '^Tf 
holds. That is to say, upon growth of Koch network 
from t generation to generation t -I- 1, the trapping time 
to first reach the trap increases by a factor of 4, which is 
consistent with Eq. (|^^ . This scaling relation is a basic 
character of the trapping process on the Koch network, 
which will be useful for deriving the formula of MTT that 
will be given in the following section. 

Notice that the order of matrix (L^)^^ is {Nt — 1) x 
{Nt — 1), where Nt increases exponentially with i, as 
shown in Eq. ([T]). Thus, for large the computation of 
the TT and the MTT from Eqs. (gOD and dH]) is pro- 
hibitively time and memory consuming, making it diffi- 
cult to obtain Tf and {T)t through direct calculation for 



large network; one can compute directly the MFPT only 
for the first several generations. However, the recursive 
construction of Koch network allows one to compute an- 
alytically the MTT to achieve a closed-form solution; the 
derivation details of which will be given in next section. 

B. Analytical solution for mean trapping time 

We now determine the average of the mean time to 
absorption, aiming to derive an exact solution. We rep- 
resent the set of nodes in Kt as Af and denote the set 
of nodes created at generation t by At. Thus we have 
A( = At U At_i. For the convenience of computation, we 
define the following quantities for m < t: 

-^m.sum ~ ^ ^ "^i (42) 
teA,„ 

and 

T* = V T* (43) 

m,sum / J I \ 

Then, we have 

rpt — I (A A) 

t,sum t — l,sum ' f,sum' V / 

Next we will explicitly determine the quantity T^^^^. To 

this end, we should first determine T^.^^^^. 

We examine the mean time to absorption for the first 
several generations of the Koch network. Obviously, for 
alH > 0, rf = 0; for t = 0, it is a trivial case. We have 
= T3 = 2. In the case of i = 1, by construction of the 
Koch network, it follows that = ^{l+Tl) + l{l + Tl), 
n = i(l + ri) + i(l+Ti), Tgi = l{l+Tl) + \{l + T^)^ 
T,i = \{l + Tl) + \{l+Tl), Ti = l{l+T^) + l{l + T^), 
and Ti ^ i(l + Tl) + i(l + Ti). Thus, 

Tlsum = E = ^4 + Ti + Ti + Ti + Ti + Ti 
ieXi 

= 12 + 2{Ti + Ti + Ti) = 12 + 2tI,^^ .(45) 
Similarly, for t = 2 case, we have 

33 

7p2 _ \ ^ rp2 _ \ ^ rp2 

^ 2, sum — 2-^ i ~ 2^ i 

tGA2 «=10 

= 48 -f 2^ {Tl + r| + r|) 

+2 {Tf + + Tl + T^ + Tl + Tl) 
= 12x4i + 22To^,,,„ + 2Tj,3„,„. (46) 

Proceeding analogously, it is not difhcult to derive that 

^t,sum = 12x4 +2 g^jjj -I- 2 T(_2,sum + ' ' ' 

+ 2*-iT*^,,,„-|-2*T*,_„ (47) 
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TABLE I; Trapping time T* for a random walker starting from node i on Kt for various t. Notice that owing to the obvious 
symmetry, nodes in a parentheses are equivalent, since they have the same trapping time. All the values are calculated 
straightforwardly from Eq. (|40|l . 



t\i 


(2,3) 


(4,5) 


(6,7,8,9) (10,11,12 





2 








1 


8 


2 


10 




2 


32 


8 


40 


2 


3 


128 


32 


160 


8 


4 


512 


128 


640 


32 


5 


2048 


512 


2560 


128 


6 


8192 


2048 


10240 


512 



34 10 42 

136 40 168 

544 160 672 

2176 640 2688 

8704 2560 10752 



TABLE IL Mean trapping time obtained by direct calculation 
from Eq. (|4l|. 



t 


Nt 


V^iVt rpt 







3 


4 


4/2 


1 


9 


60 


60/8 


2 


33 


896 


896/32 


3 


129 


13888 


13888/128 


4 


513 


219648 


219648/512 


5 


2049 


3501056 


3501056/2048 


6 


8193 


55951360 


55951360/8192 



and 

^t+l,sum = 12 X 4* + 2Tj gm,j + 2^ T(_i,sum + ' ' ' 

+ 2* g^jj^ + 2*+^ , (48) 

where 12 x 4*^"'^ and 12 x 4* are indeed the double of node 
numbers, which were generated at generations t and i+1, 
respectively. Equation (|48p minus Eq. (|47p times 8 and 
making use of the relation T^^^ — 4:Tl, one gets 

TZlsnm-^'^x^' = 2T;+^„,+8(T;^_ -12x4*-i), (49) 
which may be rewritten as 

T;+;,,,^ = 16T;^,,^-12x4*. (50) 
Using gjj^ — 44, Eq. ([50]) is solved inductively as 

T;,,_ = 5x24*-i+4*. (51) 
Substituting Eq. for t]^^^^^ into Eq. ([44]), we have 

^t*sum = ^t-l,sum + 5x2"** ^ + 4* 

= 4T/riV.m + 5x2'*"'+4*. (52) 

Considering the initial condition Tq^^^ = 4, Eq. ([5^ is 
resolved by induction to yield 

TMum = |(10x4* + 3t + 2). (53) 



Plugging the last expression into Eq. (j4T|) . we arrive at 
the accurate formula for the average of the mean time to 
absorption at the trap located at node 1 on the tth of the 
Koch network: 

1 ^' 1 t 

^ 1=2 ' 

= |x4. + i + i. ,54) 

We continue to show how to represent mean trapping 
time as a function of network order, with the aim of ob- 
taining the scaling between these two quantities. Recall- 
ing Eq. ((H), we have 4* = and t = log^iNt - 1) - 5. 
These two relations enable us to write Eq. ([5^ as 

{T)t = liNt + l log^iNt - 1) + ^ , (55) 

from which it is easy to see that for large network (i.e., 
Nt 00), the following expression holds: 

{Th^lm. (56) 

Thus, the mean trapping time grows linearly with in- 
creasing order of the network, which is consistent with 
the conclusion obtained in the preceding section. 

We have checked our analytical formulas, i.e., 
Eqs. ([54|) -([56 ll . against numerical values quoted in Ta- 
ble |TT1 For the range of 1 < i < 6, the values obtained 
from Eq. ([54]) or Eq. ([55]) completely agree with those 
numerical results on the basis of the direct calculation 
through Eq. (|¥T|l (see also Fig. [7] for comparison) . This 
agreement serves as an independent test of our theoreti- 
cal formulae. 

Notice that this linear scaling between the average 
trapping time and the network order has been previously 
obtained for three-dimensional regular lattice by using a 
method of generating function [S^l- It is also interest- 
ing to stress that this linear scaling is in contrast to the 
sub-linear scaling of mean trapping time obtained for the 
two-dimensional Apollonian network (601] and the pseud- 
ofractal scale- free web ^67|] , in spite of the fact that they 
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T 1 1 ' r 




T 1 1 ' 1 ' 1 1 r 

2 4 6 8 lU 



t 

FIG. 7: Mean first-passage time (r)t versus f on a semiloga- 
ritlimic scale. Tlie solid line corresponds to the relation be- 
tween {T)t and network order Nu (r)t « | iVt = | (2 x 4* -f 1) 
as given by Eq. (f55)) . 




FIG. 8: Second construction method of Koch network that 
highlights self-similarity. The graph after t + 1 construction 
steps, K(t -1-1), is composed of four copies of K{t) denoted 
as K1 (0 — 1, 2, 3, 4), which are connected to one another as 
above. 



have a similar topological structure [stI. IssI. IsqI. [6^. Issl. [63| 

as that of the Koch network. The reason for this dispar- 
ity is worth studying in the future. 



VII. CONCLUSION 

In this paper, on the basis of the well-known Koch frac- 
tal we have proposed a scale-free network that is called 
Koch network. We have provided a detailed exact anal- 
ysis of the topological features. We have shown that 
the Koch network displays a rich structural behavior: it 
is simultaneously scale free and small world, and has a 
high clustering coefficient. In particular, we have shown 
that the Koch network is an absolutely two-point uncor- 
related network. The especial structural characteristics 
make Koch network unique within the class of scale-free 
networks. 

The uniqueness of the topological features for Koch 
network presented here makes it potentially interesting to 
study dynamical processes occurring on the network. We 
have investigated the random walk process and a particu- 
lar trapping problem (with a trap located at a hub node) 
on the Koch network. We have obtained analytically the 
scaling efficiency exponents which are interesting to ran- 
dom walks. We have shown analytically that the mean 
first-passage time behaves linearly with the number of 
network nodes [see Eqs. (|39)) and ([56|) ]. which indicates 
that despite the presence of loops in the Koch network, 
the linear scaling of the MFPT for random walks is sim- 
ilar to that of a scale- free tree [4l|. Our analytical re- 
sult confirms a previous conclusion obtained by numer- 
ical simulations that for scale-free networks (with loops 
or not): the MFPT of random walks increases linearly 
with the network order [s^, |4l| . 
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APPENDIX A: DERIVATION OF THE AVERAGE 
PATH LENGTH 

The Koch network has a self-similar structure that al- 
lows one to calculate dt analytically. The self-similar 
structure is obvious from an equivalent network construc- 
tion method: to obtain K{t-\-l), one can make four copies 
of K{t) and join them at the hub nodes. As shown in 
Fig. [51 network K{t + 1) may be obtained by the juxta- 
position of four copies of K{t), which are labeled as K}, 
Kf, Kf, and Kf, respectively. Then we can write the 
sum Dt+i as 

A+i=4A + At, (Al) 

where At is the sum over all shortest paths whose end 
points are not in the same Kt branch. The solution of 
Eq. ((JT|l is 

t-i 

Dt = 4*-iDi+^4*-^-iA^. (A2) 

x=l 

The paths that contribute to At must all go through at 
least one of the three edge nodes (i.e., X, Y and Z in 
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Fig. [8]) at which the different Kt branches are connected. 
The analytical expression for A^, called the length of 
crossing paths, is found below. 

Denote A"''' as the sum of length for all shortest paths 
with end points in and , respectively. If K" and 
iff meet at an edge node, A"''' rules out the paths where 
either end point is that shared edge node. For example, 
each path contributing to A^'^ should not end at node 
X. If and K^ do not meet, A"''^ excludes the paths 
where either end point is any edge node. For instance, 
each path contributive to Aj ''^ should not end at node X 
or Y . Then the total sum Af is 



A, =A. 



1.2 



1,4 



A, 



2.3 



2,4 



A 



3.4 



(A3) 



By symmetry, Aj 
A(' , so that 



A^ 



A^'' and At^'^ 



Af 



and 



,2,3 



E 



X 



l-XY 



= 2{Nt - 1) St + {Nt - 1)\ (A9) 
where dxv = 1 has been used. Substituting Eqs. 
and jMI into Eq. jXH, we obtain 



Af 3Aj'^ + 3Af^ . 
In order to find A]'^ and A^'^, we define 

St = ^ dix ■ 

ieK(f) 

i^X 



(A4) 



(A5) 



Considering the self-similar network structure, we can 
easily know that at time t + 1, the quantity St+i evolves 
recursively as 



At - 12(7Vt - 1) St + 3 (A^t 
= (12< + 27) X 16*. 



If 



(AlO) 



st+i = 2st + [st + {Nt-l)] + [st + {Nt-l)] 

= 4 St + 4 X 4*. (A6) 

Using si = 12, we have 

st=4*(i + 2). (A7) 

On the other hand, by definition given above, we have 



Inserting Eq. (|A10p for Aa, into Eq. 
72, we have 



2]) and using Di 



E 
E 



{dix + djx) 



= {Nt - 1) 51 d^x + (At 

ijiX 

= 2{Nt - 1) 51 d^x 

i^X 

= 2{Nt-l)st, 



1) E 



A = IT X (2 + 7 X 4* + 6t X 4*) 



(All) 



Inserting Eq. (jAlip into Eq. one can obtain the 

(A8) analytical expression for dt as shown in Eq. (|15p . 
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